Double Counting in LDA+DMFT - The Example of NiO 
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An intrinsic issue of the LDA+DMFT approach is the so called double counting of interaction 
terms. How to choose the double-counting potential in a manner that is both physically sound 
and consistent is unknown. We have conducted an extensive study of the charge transfer system 
NiO in the LDA+DMFT framework using quantum Monte Carlo and exact diagonalization as 
impurity solvers. By explicitly treating the double-counting correction as an adjustable parameter 
we systematically investigated the effects of different choices for the double counting on the spectral 
function. Different methods for fixing the double counting can drive the result from Mott insulating 
to almost metallic. We propose a reasonable scheme for the determination of double-counting 
corrections for insulating systems. 



I. INTRODUCTION 

The combination of the density functional theory 
(DFT/LDA), a model Hamiltonian and the dynami- 
cal mean field approximation (DMFT)i, a methodology 
commonly referred to as LDA+DMFT, is to date one 
of the best approaches for the realistic description of 
strongly correlated electron systems^. While density 
functional theory does not include all the interactions 
between strongly correlated d or / electrons, it captures 
some portion of them through the Hartree and exchange- 
correlation terms. By introduction of a model Hamilto- 
nian into the calculations one tries to account for as much 
of the interactions as possible through the Coulomb in- 
teraction matrix of the impurity model. This ultimately 
leads to the problem that some contributions to the inter- 
action are included twice. This has to be explicitly com- 
pensated by adding a shift in the chemical potential of 
the correlated orbitals to the Hamiltonian, leading to the 
prominent issue of double counting. The LDA+DMFT 
Hamiltonian can be written as follows 
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Figure 1: Schematic illustration of the effect of the Coulomb 
interaction on the energy levels in a Mott-Hubbard (a) and a 
charge transfer insulator (b). Figure frorrA 



the double counting in LDA+DMFT on the spectrum 
are seldom performed. In the work presented here we at- 
tempt to shed some light on the double-counting problem 
using the example of nickel oxide (NiO). In recent years 
a number of authors applied the LDA+DMFT method 
in different flavors to this system generating a body of 
promising results- - — . 



II. NIO - A CHARGE TRANSFER SYSTEM 



E 



u, 



mm m m 



imcr im' a' 



2, (TO"' .mm' m" m" 



O VfTV (7 



creates a 



where Hlda is the LDA Hamiltonian, 
particle with spin a in a localized orbital m at site i 
and Vmra'm"m"' is the Coulomb interaction matrix be- 
tween localized orbitals. Above Hamiltonian contains the 
double-counting correction 



Hdc — Me2c E ' 



where n mj(T = c^Cma and fi^c is the double-counting 
potential. How to choose the double-counting potential 
in a manner that is physically sound and consistent is 
unknown and systematic investigations of the effects of 



Nickel oxide is a strongly correlated transition metal 
oxide that is a prototypic member of the class of charge 
transfer insulators. According to Zaanen, Sawatzky and 
Allen transition metal oxides can exhibit a behavior dif- 
ferent to the classic Mott-Hubbard picture^. In a Mott- 
Hubbard insulator the charge gap opens through splitting 
of the d band by the Hubbard U. In the charge-transfer 
system the gap typically opens between hybridized lig- 
and p and transition metal d states and the upper Hub- 
bard band corresponding to the d states of the transition 
metal. Thus, it is not only the Hubbard U, but also the 
so called charge transfer energy A = |e<j — e p \ that de- 
termines the size of the gap. In the scheme by Zaanen, 
Sawatzky and Allen materials can be classified by their 
respective values of U and AiS. For A > U the system 
is a Mott-Hubbard insulator, whereas for A < U it be- 
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longs to the charge transfer class. In general, systems 
with completely filled d(e g ) and partially filled d(t 2g ) 
shells, like titanates, vanadates and some ruthenates be- 
long to the Mott-Hubbard class. Prominent examples of 
charge transfer insulators are NiO, MnO, manganites and 
cuprates. In these systems the e g shell is partially filled 
and the t 2g shell is fully occupied. 

The density of states and the band structure of NiO as 
obtained by LDA (using the PAW±i based VASP code^) 
are shown in Fig.©. The band structure shows five Ni 
3d bands in the energy window — 2.5eV to +1.5eV cross- 
ing the Fermi energy and three separated O 2p bands 
below, extending down to — 8eV. These bands contain 14 
electrons in total, 6 occupy the oxygen p bands and the 
remaining 8 the Ni d bands. In contrast to the LDA pre- 
diction NiO is not a metal, on the contrary, experiments 
revealed a charge gap of about 4eV— . Additionally, it 
exhibits antiferromagnetic order below the Neel temper- 
ature of TV = 525K. Our computations were carried 
out in the paramagnetic phase, which is not problematic, 
since the gap opened by electronic correlations does not 
depend on whether the system is magnetically ordered. 
It has been shown in angle-resolved photoemission ex- 
periments (ARPES), that passing the Neel temperature 
does not qualitatively alter the valence band spectrum^. 



III. METHODOLOGY AND RESULTS 

The model that has to be used for a simulation of 
NiO is the five band Hubbard model which describes 
the correlated 3d states of Ni. We have calculated the 
model parameters of such a model in an ab initio fashion. 
The local orbitals are represented by Wannier functions, 
which have been shown recently 15-17 to be a very good 
choice for a basis set, because they form a complete ba- 



sis of the Hilbert space spanned by Bloch functions and 
are reminiscent of localized atomic orbitals. Our calcula- 
tions involved two different flavours of the LDA+DMFT 
framework: One uses a projection of Bloch states on lo- 
cal orbitals represented by Wannier function a 17 ' 18 and a 
Quantum Monte Carlo (QMC) solver—, while the other 
employs the Linear Order Muffin-Tin Orbital method 
(LMTO) ?0 and a finite-temperature exact diagonaliza- 
tion (ED) solver^ 2 .. 

The effective Wannier Hamiltonian includes the five 3d 
bands of nickel as the correlated subspace and the three 
2p bands of oxygen as the uncorrelated part. The in- 
clusion of the p bands is physically motivated since in a 
charge transfer compound the oxygen bands play an im- 
portant role in the physics of the system, as was pointed 
out above. A computation taking into account only the 
Ni d states is capable of reproducing the insulating be- 
havior and the size of the gap as shown by Ren et alA 
Additionally, the double counting is reduced to a trivial 
shift in calculations that contain only the Ni d bands, 
since the full Wannier Hamiltonian belongs to the cor- 
related subspace. The double counting can thus be ab- 
sorbed into the total chemical potential. However, the 
physics of the charge transfer insulator cannot be cap- 
tured without taking into account the ligand p states. 

Our calculations were performed at inverse tempera- 
ture (3 = 5eV _1 , which corresponds to 2321K, using 
up to 80 time slices and on the order of ~ 10 6 Monte 
Carlo sweeps in the QMC. In the ED fraction of calcu- 
lations we used a ten site cluster (5 impurity levels and 
5 bath levels). The temperature used may appear high, 
yet it is low enough to give a qualitatively correct de- 
scription of the physics of the material. Computations 
at lower temperatures pose no fundamental problem, the 
amount of LDA+DMFT calculations performed for this 
study would have made them too expensive though. We 
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Figure 3: Spectral functions at f3 = 5eV 1 for different values of the double counting fide obtained with LDA+DMFT (QMC). 



have used a Coulomb interaction matrix corresponding 
to the parameter values U = 8eV and J = lcV. 

The double-counting potential fide defined above is 
found to have profound impact on the density of states 
Ni(u>) = —^Gi(uj) shown in Fig. ^ and the k-resolved 
spectral function 

Ai(k,uj) = — ilm (w + /i — £j(k) — Ej(oj)) 1 shown 
along the line T — X in the Brillouin zone in Fig. ((3]). The 
spectral functions were obtained by the maximum en- 
tropy method^ from imaginary time Green functions. 
The double-counting potential has been treated here as 
an adjustable parameter and has been varied between 
21eV and 26eV. These values already contain the intrin- 
sic shift due to the energy of the particle-hole symmetry 
in the Hirsch-Fye QMC method that amounts to 34eV 
with our values of U and J. The energy of the particle- 
hole symmetry is obtained from Eq.(fT]) with n° = h. 

The most prominent effects of the double counting 
on the spectral properties are the shift of the oxygen 
p bands with respect to the nickel d bands, as well as 
the variation in gap size. Plainly speaking, the double- 
counting correction allows for a tuning of the spectral 
properties from a large gap Mott-Hubbard insulator to 
a metal. The regime of the charge transfer insulator, 
the expected physical state of NiO, lies somewhere in be- 
tween. The experimental spectrum, obtained by x-ray- 
photoemission (PES) and bremsstrahlung-isochromat- 
spectroscopy (BIS) showing both occupied and unoccu- 
pied parts, was obtained by e.g. Sawatzky and Alleni^. 
The spectrum recorded at 120eV is predominantly of Ni 
3d character, while the 66eV spectrum contains a strong 
contribution of O 2p at about — 4eV— Additionally, 
the detailed decomposition of the spectra showed contri- 
butions of both O 2p and Ni 3d at the top of the valence 
bandit. The calculated LDA+DMFT(QMC) spectral 
function shown in Fig. ^ show basically the two differ- 
ent physical situations of a Mott-Hubbard Fig. (3(a) I and 
a charge-transfer insulator Fig.( |3(b)"j ) mentioned above. 
Both spectral functions were obtained for NiO, by vary- 
ing the double-counting correction. The characteristic 



feature of a charge-transfer system, the strongly hy- 
bridized ligand p and transition metal d character of the 
low-energy charge excitations^^, is only present in the 
spectrum in Fig.( |3(b)"j ). The spectrum in Fig. (3(a)) is 
missing this feature almost completely and shows Mott- 
Hubbard behaviour. This difference underscores the im- 
portance of the proper choice for the double-counting cor- 
rection. 

Let us now turn to the k-resolved spectral func- 
tions shown in Fig.@ and compare them with ARPES 
data-£i2S. The uppermost band in Figs. ( 4(a)] |4(b)| at 



2eV above the Fermi level is a Ni e g band, while the 
other bands can be identified with the ones obtained 
by ARPES. The two lowest lying bands correspond to 
oxygen p states, the bands above are formed by Ni d 
states. The characteristic features seen in ARPES, like 
the broadening of the oxygen bands around the midpoint 
of the r — X line, are clearly present. The quantita- 
tive features, especially the relative band energies can 
strongly differ, depending on the double counting cho- 
sen. The bands in Fig. (4(a) I (fide — 21eV) show a clear 



separation between the oxygen and the nickel part at the 
T-point as well as the X-point. At the increased value 
of the double counting fide = 25eV the oxygen bands are 
shifted towards the Fermi level, coming to overlap with 
the Ni d bands at the T point as in the ARPES data. 
A detailed comparison of the calculated bandstructures 
with experiments shows that the bands calculated with 
fide = 25eV agree very well with the experimental data. 
These calculations reproduce the flat bands at — 4eV and 
another at about — 2eV becomes more prominently visi- 
ble at fide — 25eV, while it is very faint at fide = 21eV. 
The dispersive bands in the region — 4eV to — 8eV also 
agree very well with experiment. Our calculations at this 
value of fide yield very similar results as those obtained 
by Kunes et al£. Calculations with other values of the 
double counting can strongly differ from the experimen- 
tal data, as shown by the example of fide = 21eV in 
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Figure 4: k-resolved spectral functions A(k,u>) along the line V — X in the Brillouin zone for different values of the double 
counting fi dc obtained using LDA+DMFT (QMC). 



The dimension of the problem of the double counting 
becomes apparent if the parameter space of the overall 
chemical potential fi and the double-counting potential 
fide versus the total particle number in the system N is 
examined. The result is shown in Fig. ([5]) with the particle 
number color coded. The picture shows that in princi- 
ple any combination of fi and fide that yields a point in 
the green plateau, corresponding to the desired particle 
number TV = 14 a priori describes the system equivalently 
good. The problem that arises here is that convention- 
ally fixing the total chemical potential fi in the middle 
of the gap still leaves one the freedom of choosing differ- 
ent values for fi dc - An additional condition is required to 
completely determine the systems position in the (fi, fide) 
parameter space and thus in the end its spectral prop- 
erties. As we have argued above this choice is of crucial 
relevance for the results of the LDA+DMFT simulation 
and not just an unimportant technicality. 

Since other, related approaches, like the LDA+U 
method, also include a double counting the problem is 
not new. Over the years different analytic methods to fix 
fide have been devised. Two prominent examples are the 
around mean-field (AMF) 27 approximation and the fully 
localized or atomic limit (FLL)2&. The AMF is based 
on the conjecture that LDA corresponds to a mean-field 
solution of the many-body problem, as was argued by 
Anisimov et al£l. The resulting double-counting poten- 
tial is 



AMF _ rj o 



J mm')^ ) (1) 



where n 



2(21+ 
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is the average occupancy. 



We use the global average and not the spin dependent 
version proposed in Refi^, since we were performing 
paramagnetic calculations in which both spin compo- 
nents are equally occupied. One assumes all orbitals 
belonging to a certain value of the angular momentum 
I to be equally occupied and subtracts a corresponding 
mean-field energy. This is, however, incorrect, since LDA 
contains the crystal field splitting explicitly and will in 
general not produce equally occupied orbitals even for 
weakly correlated systems. The result for the case of NiO 



using self-consistent occupancies from the DMFT loop is 
shown in Fig.© labeled (SC)AMF. The value obtained 
with the formal occupancies given above ((F) AMF) lies 
outside of the considered part of the parameter space 
at 20.4eV. In both cases the solution corresponds in ou r 
case to a Mott- Hubbard insulator as shown in Fig. (3(a) I. 
The AMF functional is known to produce unsatisfactory 
results for strongly correlated systems, which led to the 
development of another method, the so called FLL. 

The FLL functional takes the converse approach to 
the AMF and begins with the atomic limit. It has been 
shown, that this new potential can be written as a cor- 
rection of the AMF solution ([]} in the following for 
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H F dc LL =H A d r + {U-J){n a -^). 

This addition to the AMF potential has the effect of a 
shift of the centroid of the level depending on its occupa- 
tion. An empty level is raised in energy by ^(U — J) and 
the converse happens to a fully occupied level. The form 
of the functional is based on the property of the exact 
density functional that the one electron potential should 
jump discontinuously at integer electron number—, which 
is not fulfilled in LDA or GGA. Ultimately the FLL leads 
to a stronger trend towards integer occupancies and lo- 
calization. The result of the FLL, as shown in Fig.©, 
constitutes a substantial improvement over AMF, yet still 
produces too low values. The general problem with an- 
alytic expressions like the ones presented is that their 
scope is limited to certain classes of systems that ful- 
fill the assumptions made in the derivation process. The 
AMF for example might give good results for weakly cor- 
related systems, but it certainly fails for the strongly cor- 
related ones. The FLL improves the situation for insula- 
tors, but it is still based on ad-hoc assumptions. Addi- 
tionally a certain degree of ambiguity is inherent, since 
one can compute the corrections using the formal occu- 
pancies given above, occupancies obtained from LDA or 
from the self-consistent DMFT loop. Other analytical 
formulas for the double-counting correction have been 
proposed for the case of NiO, see e.g. the work by Ko- 
rotin et al^ and Kunes et al£. Despite giving reasonable 
resulting spectral functions analytical approaches to the 
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Figure 5: Surface created by different combinations of the chemical potential /i and the double-counting potential /j,d c plotted 
versus the particle number N obtained with LDA+DMFT (QMC). The particle number has been color coded: the green plateau 
corresponds to a particle number very close to the desired value of N = 14, values below are encoded in blue, values above in 
red. Additionally the results produced by different methods to fix the double counting are indicated. For the AMF and FLL 
functionals SC or F in parentheses indicates, that the self-consistent occupancies from the DMFT or the formal occupancies 
have been used respectively. For further details we refer to the text. 




double counting are not optimal. 

The obvious problems with analytical formulas make 
conceptually different approaches worth exploring. It 
would certainly be an improvement if the double count- 
ing could be found self-consistently along with the chem- 
ical potential in the DMFT self-consistency loop. Since 
the double counting correction is intrinsically an impu- 
rity quantity and not a global quantity (like the chemical 
potential fi) it would be most desirable to use intrinsic 
quantities of the impurity like the impurity self-energy or 
the impurity Green function to fix it. One possible ansatz 
using the impurity self-energy S^^, is to constraint the 
high energy tails in the real part of the self-energy to sum 
up to zero 

ReTr(E3^ JV ))=0. 

Here, a; at is the highest Matsubara frequency included in 
the computation. Physically this amounts to the require- 
ment that the shift in the centroid of the impurity orbitals 
contains no static component. The resulting correction is 
fide ~ 21.3eV and thus very close to the (SC)AMF value 
shown in Fig.©. The result produced is thus reason- 
able in principle in the sense that it produces an insulat- 
ing solution. However, the resulting spectrum resembles 
a Mott-Hubbard system. Double counting corrections 
based on the self-energy have been applied successfully 
to metallic systems, see e.g<2£. 

Another very promising approach, which is in princi- 
ple based on the Friedel sum rule^i, is to constraint the 
total charge in the impurity. This approach requires that 
the electronic charge computed from the local noninter- 
acting Green function and the one computed from the 



interacting impurity Green function are identical 18 

TrG3(/3)=TrG° m ^(/3). (2) 

Alternatively one can also use the Weiss field Q mm ' in- 
stead of the local noninteracting Green function in above 
equation. Both versions of the method give very similar 
results and work very well in metallic systems^ 8 -, since 
in a metal the total particle number of the system N 
and of the impurity rii mp are both very sensitive to small 
variations in fj, and fide- As NiO has a quite large gap 
the charge does almost not vary with neither the chem- 
ical nor the double-counting potential in the gap. The 
constraint of fixed particle number can thus be fulfilled 
to a very good approximation in the whole gap region, 
the criterion j2]) essentially breaks down. Since the gap 
in NiO is large this method fails and drives the system 
towards a metallic state at double counting fj, ( ic ~ 26.5eV 
indicated by the arrow pointing out of Fig.©. 

Since the double-counting corrections that we have ex- 
plored either fail to reproduce the physics of NiO or are 
based on analytic arguments that do not exactly apply to 
the system a different, sound way fixing the value of the 
double counting for insulating systems is needed. Since 
the double-counting potential effectively acts like an im- 
purity chemical potential we propose to find the value 
at which it lies in the middle of the gap of the impurity 
spectral function where the occupation of the impurity 
is about riimp sa 8 particles. This part of the calcula- 
tions was done using the exact diagonalization impurity 
solver (see above), which is much faster and uses the 
full Coulomb interaction matrix including spin-flip and 
pair-hopping terms. Additionally it does not suffer from 
statistical errors and directly provides data on the real 
axis. We used a 10 site cluster with 5 impurity levels plus 
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5 bath levels and fit the bath Green function via the level 
energies and hopping parameters^. An explicit scan of 
the parameter space revealed that the proper value for 
the double-counting correction is /idc ~ 25eV, basically 
the same value found above by inspection and compari- 
son of spectral features to experimental data. It is indi- 
cated as INS in Fig.©. The corresponding lattice and 
cluster spectral functions are shown in Fig.©. The pro- 
posed criterion thus produces a double-counting correc- 
tion that reproduces the spectral features of the valence 
in accord with photoemmission measurements and does 
not contain ad-hoc assumptions about the system. 

IV. SUMMARY AND CONCLUSIONS 

In summary, our study has shown that the double- 
counting correction in the LDA+DMFT formalism has 
to be very carefully assessed when performing calcu- 
lations with a correlated and uncorrelated part in the 
Hamiltonian. We have examined the influence of the 



double-counting potential on the spectral properties us- 
ing the example of NiO. Different tracks in the search 
for a sound double counting were explored. A well de- 
fined analytical expression for the double-counting po- 
tential fide probably cannot be formulated in the context 
of LDA+DMFT. Thus, one has to resort to numerical 
criteria to fix the value of the double-counting correc- 
tion. For metals the self consistency criterion based on 
the charge Eq.© works very reliably. It is, however, not 
applicable to insulating systems. In such a case we pro- 
posed to fix the value of the double-counting potential by 
setting it in the middle of the gap of the impurity spec- 
tral function. This criterion led to spectral properties 
in good agreement with experiments. Thus, one has to 
resort to self-consistent numerical approaches to fix the 
double-counting correction properly. Further work, espe- 
cially the examination of other systems will show if the 
proposed methodology can be reliably applied to predict 
the electronic structure of correlated electron systems by 
LDA+DMFT calculations. 
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